Isoscaling behavior in the Isospin dependent Quantum Molecular Dynamics model 
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The isoscaling behavior is investigated in the frame of Isospin dependent Quantum Molecular 
Dynamics (IQMD) models. The isotopic yields ratio Y 2 /Yi for reactions 4S Ca+ 4S Ca and 40 Ca+ 40 Ca 
at different entrance channels are simulated and presented, the relationship between the isoscaling 
parameter and the entrance channel is analyzed, the results show that a and /3 reduce with the 
rise of incident energies and increase with the impact parameter b, which can be attributed to the 
temperature varying of the pre-fragments in different entrance channels. The relation of a and 
symmetry-term coefficient C sym reveals that the chemical potential difference Afj, is sensitive to the 
symmetry-term coefficient C sym . 



PACS numbers: 25.70. Pq, 24.10 Nz 

Experimental investigations of the properties of isospin 
asymmetric nuclear matter are receiving increased at- 
tention due to the possibility of creating nuclear mate- 
rial with appreciable isospin asymmetry using neutron- 
rich beams. This is expected to enable one to extrap- 
olate the present understanding regarding the effective 
nucleon-nucleon interaction to the unknown domain of 
large isospin asymmetry. The isospin asymmetry of re- 
action parts plays an important role in the evolution of 
nucleon exchange processes which form the basis of many 
nuclear transport models 0, 0, IS 0, IE IE ■ It is also an 
effective indicator of the degree of chemical equilibrium 
reached in heavy ion reactions. 

To minimize undesirable complications stemming from 
the sequential decays of primary unstable fragments, it 
has been proposed that isospin effects can be best stud- 
ied by comparing the same observable in two similar reac- 
tions that differ mainly in isospin asymmetry [H l§l liSl fiof . 
If two reactions, 1 and 2, have the same temperature but 
different isospin asymmetry, for example, the ratio of a 
specific isotope yields with neutron and proton number N 
and Z obtained from system 2 and system 1 have been ob- 
served to exhibit isoscaling, i.e., exponential dependence 
of the form 



R 2 i(N, Z) = 



Y 2 (N, Z) 
Yi(N,Z) 



Cexp(aN + (3Z). (1) 



where a and (3 are the scaling parameters and C is an 
overall normalization constant, with the convention of 



* Supported by the Major State Basic Research Development Pro- 
gram under Contract No G2000774004, the National Natural Sci- 
ence Foundation of China (NNSFC) under Grant No 10328259 and 
10135030, and the Chinese Academy of Sciences Grant for the Na- 
tional Distinguished Young Scholars of NNSFC. 
t Email: wdtian@sinap.ac.cn, Tel:021-59556905 
i Corresponding author: ygma@sinr.ac.cn 



that the neutron and proton composition of reaction 2 is 
more neutron-rich than that of reaction 1. The system- 
atization of the experimental data in form (1) has been 
conformed to variety of reactions over wide energy range 
and reaction systems [jj 0, . This scaling law occurs 
naturally within various reaction processes, such as the 
evaporation, deep inelastic scattering fission [T^.[T^|. 
and multifragmentation [jj, 0, . In the framework of 
different models |9| , it has been revealed that the isotope 
yields ratio are not affected very much by the sequen- 
tial decays following the fast stage of the reaction, so 
R,2i(N, Z) deduced from the detected isotope fragments 
can reflect the isotope yields ratio of the primary frag- 
ments. The isoscaling parameters a and (3 can be used 
to study the isospin dependent properties in nuclear col- 
lisions, and this method allows one to study the early 
reaction stage of the decay fragmenting system, to un- 
derstand the dependence of measured isotope distribu- 
tion on properties of hot emitter, the nuclear asymme- 
try term in nuclear matter equation of state (EOS), and 
allows the direct comparison between experimental mea- 
surement with theoretical simulation. 

To explore the isospin properties and the fragmenta- 
tion in dynamical nuclear collisions, we study the isoscal- 
ing phenomenon by the IQMD model IToL flTL fl^| , which 
is based on the general QMD [Tfjl I20L I2llf22fl model to 
include explicitly isospin-degrees of freedom. The QMD 
model is classical in essence because the time evolution of 
the system is determined by classical canonical equation 
of motion, however, many important quantum features 
are included in this prescription. It is well known that 
the dynamics in heavy ion collision at intermediate ener- 
gies is mainly governed by three components: the mean 
fields, two-body collision and Pauli blocking, therefore, 
for isospin-dependent reaction dynamics model it is im- 
portant for these three components to include isospin de- 
grees of freedom. In additional, in initialization of pro- 
jectile and target nuclei, the samples of neutrons nd pro- 
tons in phase space are also treated separately since there 
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exists a large difference between neutron and proton den- 
sity distributions for nuclei far from the /3-stability line 

HlUzl- 

In the IQMD model, the nuclear mean field is given by 



U(p, t z ) = U dd + U Yuk + U Coul + U Sym + U 



MDI 



(2) 



where U dd is the density-dependent (Skyrme) poten- 
tial, U Yuk the Yukawa (surface) potential, JJ Cow[ the 
Coulomb energy, U Sym the symmetry energy term and 
jjMDi ^. ne m0 nientum dependent interaction, but in our 
present calculation, JJ MDI was not included. The U dd 
can be written as 



U dd = a{-^) + b(£-)i 
Po Po 



(3) 



with po — 0.16/m 3 (the normal nuclear matter density), 
a = -356 MeV, b — 303MeV, and 7 = 1.17 (corresponding 
to the Soft EOS). 



TjYuk "V 

*F3 



-exp(Lm 2 ) 



x [exp(mrij)er f c(VLm — /y/A~L) 
—exp(mrij)erfc(VLm + rij/VAL)] (4) 



Neutron Number N 
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Atomic Number Z 



jjcoui = e ^ _ (1 + _ T . <)er/c(ry/> /4L) (5) 

7/Sj/m - Csym V T T ■ - CXV\ ^ n ~ r ^ 2 \ (f\) 

with = -0.0024 GeV, m = 0.83, and the L is the so- 
called Gaussian wave-packet width (here L = 2.0 /m 2 ). 
The relative distance T\j — \v\ — The ri Z is the z— th 
component of the isospin degree of freedom for the i— th 
nucleon, which is equal to 1 and -1 for proton and neu- 
tron, respectively. C sym is symmetry energy strength. 
More details on QMD and IQMD models are discussed 
in Ref [H [H M HJ . 

We first perform the reaction simulations for 40 Ca + 
40 Ca and 48 Ca + 48 Ca collisions with impact parameter 
b = at several different incident energies E/A = 25, 
35, 50 and 70 MeV, the dynamical process was simulated 
until t = 200 fm/c. In central collisions, many fragments 
are formed, which can be basically divided into projectile- 
like, target-like and neck region emitted products. In the 
neck region of central collision, it is more reliable that the 
thermal and chemical equilibrium could be achieved, so 
we select the products with the parallel velocity condition 
of \v///vi n \ < 0.2, where v// means the parallel velocity 
of the products, and Vi n is the incident velocity, by which 
determined the fragments emitted from the neck region 
parts. In the following discussion all simulation results 
are under such selection. 



FIG. 1: The fragment yield ratio of the IQMD simulations for 
central collisions 48 Ca + 48 Ca and 40 Ca + 40 Ca at 35MeV/A, 
a and /3 are printed inside the figure. 



Fig. 1 shows the yields ratio Y N , z ( 48 Ca + 48 
Cd)jY^^z^ a Ca + 40 Ca) of these two reaction at the 
incident energy E/A=35MeV, clear isoscaling relations 
are observed, the error bar inside Fig. 1 and following's 
plots are statistical error only.The isoscaling parameters 
a and [3 can be extracted from formula (1), which are 
printed inside the figure with a — 0.236, and — -0.326. 
While extracting a and /3, we use R21 (N) — C'exp(aN) 
and Rix(Z) = C" exp({3Z), respectively. Fig. 2 plots 
the variation of a and (3 with the incident energies for 
central collisions 40 C a + 40 Ca and 48 Co + 48 Ca , one 
can find that the isoscaling parameters a decreases with 
the increasing of the incident energy as well as \/3\. In 
the grand-canonical approximation p|, , a = Kp, n /T 
and (3 = Ap, p /T, provided a common temperature T for 
both systems exists 0iA*n an d are neutron and pro- 
ton chemical potential of the system, determined by the 
system composite, the density p, the temperature T and 
its overall N/Z ratio [ITl Il5|. In two similar reactions 
like 40 Ca+ 40 Ca and 48 Ca + 48 Ca, Ap n and Ap, p should 
not show any significant changes with the temperature 
for light charged particles as the model calculations [lT| . 
therefor the reduction of a and \j3\ with the increasing of 
incident energies can be mainly attributed to the rise of 
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FIG. 2: The variance of isoscaling parametes a and /3 with 
the incident energies for central collisions 48 Ca+ 48 Ca and 
40 Ca+ 40 Ca 



the system temperature. 



FIG. 4: Relation between (Z/A) 2 and a for three differ- 
ent isospin systems 40 Ca+ 40 Ca, 48 Ca+ 48 Ca and 60 Ca+ 60 Ca 
with different symmetry potential coefficient C Bym , and 
the insert figure is the trend of a and C sy m of two 
yields ratios, the square solid point represents the re- 
sult from ratio of Y( 48 Ca+ 4S Ca)/Y( 40 Ca+ 40 Ca), and the 
round solid point is extracted from the yields ratio of 
Y^Ca+^C^/Y^Ca+^Ca). 
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neutron and proton chemical potential difference A/x„ 
and A/ip does not change with the varying of the tem- 
perature very much ^J, the reduction of isoscaling pa- 
rameters a and \/3\ can be caused by the temperature 
difference. 

Furthermore, to explore the aspect of equilibrium in 
fragment emission in IQMD simulation, we study the 
relationship between the isoscaling and the fragment 
isospin asymmetry in IQMD simulation. The following 
relation (7) was found relating the isoscaling parameter 
a, the (Z/A) 2 of fragments, and the symmetry energy 
coefficient C(Z) by the form [l2| 



FIG. 3: Relation of isoscaling parameters a and /3 with the 
impact parameter b, at incident energy 35MeV/A for colli- 
sions 48 Ca+ 48 Ca and 40 Ca+ 40 Ca, the dash lines are plotted 
to guide the eyes. 

In Fig. 3 the relation of isoscaling parameters a and 
|/3| with the impact parameter b is shown for collision 
48 Ca+ 48 Ca and 40 Ca+ 40 Ca, a and \/3\ values raise with 
the increasing of b. As we know, the reaction become 
less violent from central to peripheral collision, then sys- 
tem becomes less excited, the temperature T of the sys- 
tem generally be defined as T — ^Jk(Eo) / (Aq), where Eq 
is the excitation energy of the prefragment, and Aq the 
mass of the decaying prefragment, k denotes the inverse 
level density parameter, the initial temperature T de- 
termined by the excitation of the decaying prefragment. 
The peripheral collision have less excitation energy com- 
paring with the central collision, which means low tem- 
perature the emitting source can be reached. The free 



a = 4C s . tfm [(Z 1 /A 1 ) 2 - (Z 2 /A 2 ) 2 ]/T (7) 

where C syrn is the symmetry energy coefficient and T 
is the system temperature, 1 and 2 denote the neutron- 
deficient and neutron-rich reaction system, respectively. 
In IQMD simulation, three reactions of 40 Ca+ 40 Ca, 
48 Ca+ 48 Ca and 60 Ca+ 60 Ca are simulated with differ- 
ent symmetry-term potential coefficient C sym , at the in- 
cident energy E/A=35MeV, the relation between the 
(Z/A) 2 , symmetry potential coefficient C sym and the 
isoacaling parameter a from three simulations are plot- 
ted in Fig. 4, where Z/A is calculated from the emit- 
ted fragments, averaged over the fragment 5 < Z < 
8. One can find that linear relation is generally kept 
between a and (Z/A) 2 for different C sym , this linear 
function can prove the chemical potential equilibration 
in reactions 40 Ca+ 40 Ca, 48 Ca+ 48 Ca, and 60 Ca+ 60 Ca 
though different symmetry potential coefficients C sym 
are used in the simulations, since equation (7) is satis- 
fied only in thermal and chemical potential equilibrium. 
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The insert figure shows that the isoscaling parameter a 
increases with the increasing of symmetry-term coeffi- 
cient C sym , in both Y( 48 Ca+ 48 Ca)/Y( 40 Ca+ 40 Ca) and 
Y( 60 Ca+ 60 Ca)/Y( 40 Ca+ 40 Ca). Up to now no evidence 
show asymmetry-term coefficient C sym have impact on 
the system temperature, since the temperature of the 
system, as we mentioned above, mainly determined by 
the incident energy, also means the excitation energy of 
the reaction. Other models like isospin dependent lattice 
gas model (LGM) 0] and statistical multifragmentation 
(SMM) model (Tj predicted that the chemical potential 
differences between two similar reaction A/i„ varies little 
in different temperature. But the chemical potential dif- 
ference A^t n changes a lot with different symmetry-term 



coefficients C sym , the increasing of a mainly comes from 
the contribution of chemical potential difference A/i n . 

In summary, with help of the IQMD simulation, the 
isotope yield ratio in the multi-fragmentation process 
in intermediate energies heavy ion collision, shows the 
isoscaling behavior in different entrance channels. The 
isoscaling parameters a and [3 drop with the incident en- 
ergies and raise with the impact parameters, which can 
be basically attributed to the different temperature of the 
emitting source. The isoscaling law also presents a linear 
relation between a and (Z/A) 2 , but the a varies with 
symmetry-term coefficient C sym which indicates a can 
be served as a sensitive probe to the symmetry potential 
as LGM shows 15). 
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